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Abstract 

We study the hadronic production of B s and B* mesons within the fixed-flavor-number scheme, 
in which the dominant gluon-gluon fusion mechanism is dealt with by using the complete 
approach. Main theoretical uncertainties for B s and B* production at TEVATRON and LHC 
are presented. It is found that when m s increases by steps of 0.1 GeV, the integrated cross 
section of B^ decreases by 80% — 100%, when increases by steps of 0.1 GeV, it changes to be 
~ 10%. While the uncertainties caused by the parton distribution function and the factorization 
scale varies within the region of g to |. Considering possible kinematic cut on the transverse 
momentum and the rapidity cut for the detectors at TEVATRON and LHC, we also make 
estimations on the B s and B* production with various kinematic cuts. 
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Since Run II at the TEVATRON Col 



have successfully collected B s data 



ider started in 2001, the CDF and DO experiments 



One can use B s meson to study those 



interesting topics as QCD model building, physics beyond the Standard Model, the electro 



b, y, 0, 3. 



weak symmetry breaking mechanism, charge-parity (CP) violation and etc. 
Taking into account the prospects of B s production at Fermilab TEVATRON and at the 
newly running CERN LHC, the future numerous data require more accurate theoretical 
predictions, especially on its hadronic production. 

According to the QCD factorization formula, the hadronic production of B s and B* can 
be written as 



where v^S* stands for the total collision energy of the incoming hadrons, F % Hi P (xx, n 2 F ) 
and p 2 (x2, [ip) are the parton distribution functions (PDFs) of incoming hadrons Hi 
(momentum Pi) and H2 (momentum P2) for parton i (with momentum fraction x\) and 
parton j (with momentum fraction X2) respectively. Q 2 is the "characteristic energy scale 
of the subprocess squared" and p F stands for the factorization scale for the PDF and the 
hard subprocess. A detailed discussion on the choice of Q 2 and pf can be found in Ref.Q], 
here for simplicity, we shall take Q 2 = p 2 F for the present perturbative QCD calculation. 
d<7..^ B (*) x stands for the differential cross-section of the relevant hard subprocess, in which 
s = X1X2S is the c.m.s. energy of the subprocess and Pt is the transverse momentum of 

Within the fixed-flavor-number (FFN) scheme lo| . where only light quark/antiquark and 
gluon should be considered in the initial state of the hard scattering subprocess, it can be 
found that B^ hadronic production are dominated by the gluon-gluon fusion mechanism, 
which is through the sub-process g + g — > B^ + b + s and is of order a A s . In addition to the 
gluon-gluon fusion mechanism, there are several different mechanisms for the production, 
such as that via the quark-antiquark annihilation subprocess qq — > B^ + b + s and etc.. 
However, it can be found that the contributions to the production from quark-antiquark 
annihilation are much smaller (only about 1%) than those from gluon-gluon fusion, which 
is due to the fact that the 'luminosity' of gluons is much higher than that of quarks in pp 
collisions (LHC) and in pp collisions (TEVATRON), and there is a suppression factor due to 




FIG. 1: Typical Feynman diagrams for the hard subprocess g{Ki)+g(K2) — > Bi*\p)+b(qb)+s(q. 



he virtual gluon propagator in the annihilation, which is similar to the case of B c production 
Hence in the present letter, we shall concentrate our attention on the gluon-gluon 



fusion mechanism. And to be useful experimentally, we shall discuss the main uncertainties 
in estimating the hadronic production of B^ , which includes the choices of the factorization 
energy scale //p, the various versions of parton distribution functions (PDFs), the values of 
the bound state parameters and etc.. 

As for the dominant gluon-gluon fusion mechanism, its hard subprocess g + g — > B^ + 
b + s includes 36 Feynman diagrams, whose typical ones are plotted in FIG HI To derive 
the analytical squared amplitude for this subprocess is a tedious task, since it contains 
non-Abelian gluons and massive fermions. Fortunately however, very recently, a generator 
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15| for hadronic production of B c meson has been available, where to 
deal with the subprocess g + g — > B c + b + c, the so called 'helicity amplitude approach 
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171 ] 1 has been adopted to derive analytic expressions at the amplitude level and then 



do the numerical calculation just at the amplitude level. Here we adopt the same method of 
to deal with the present B^ production, which can be obtained by suitable 



changing the c-quark lines defined in Ref. 13] to the present s-quark lines. To short the paper, 
we shall only present the main idea on how to deal with those 36 Feynman diagrams based 



on the 'helicity amplitude approach', and the interesting reader may consult Ref.[l3j for 
detailed calculation technology. The main idea is to convert the problem into an equivalent 



1 It should be noted that Ref.|16lj presented only the formulae for massless spinor lines, so proper changes 
as have been done in Ref.[l3| should be made so as to deal with the massive spinor lines. 
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'massless' one that is well solved in literature, i.e. by transforming the massive quark lines 
to be massless ones, and then to apply the symmetries as much as possible. To extend 
the symmetries for the amplitude corresponding to 36 Feynman diagrams, we first focus on 
the numerator of the amplitude related to typical fermion lines, and neither consider the 
color factors nor distinguish the flavor of the fermion lines at the moment. Then, because 
of Feynman diagram symmetries, these diagrams can be grouped into a few typical ones 
according to the different type of fermion lines. And then we implement proper factors for 
the fermion lines: color factors, suitable denominator and spinors and etc., so as to obtain 
an exact and full typical fermion line that appears in Feynman diagrams. When all kinds 
of typical fermion line factors, factors for external lines of gluons and gluon propagators are 
'assembled', then the full term, corresponding to the Feynman diagram of the amplitude, 
is achieved. Next, to do the phase space integration, we first use RAMBOS 18| routine to 



generate the requested phase space points and then use VEGAS [19J program to perform 
the integrations. 

Based on the above calculation technology, we present the numerical results. As for the 
present LO estimation, fs s appears in the amplitude as a linear factor, so the production 
cross sections are proportional to it squared. Therefore, the uncertainties in the production 
from fp can be figured out straightforwardly, so throughout the paper, we take fs s = 0.209 
GeV[20|. And because the spin splitting effects are ignored here, so there is no difference for 
the decay constant between the spin stats ^Sq] and [ 3 Si]. Further more, we shall study the 
uncertainties in 'a factorization way' throughout the paper, i.e., all of the parameters vary 
independently in their reasonable regions. For instance, when focussing on the uncertainties 
from the constitute s-quark mass m s , we let it be a basic 'input' parameter varying in a 
possible range 

OAGeV <m s < OJGeV, (2) 

with all the other factors, including the f^-meson mass, the decay constant fs 3 and etc. 
being fixed. 

In TAB]U we show the uncertainties from m s , where the other factors are fixed precisely 
as: nib — 4.9 GeV, the PDFs are taking as CTEQ6L 23j; the strong coupling a s is in LO 
and the factorization energy scale is taken to be fip = p T + m% . Note that for the mass 



of B s , the experimental result is m Bg = 5.3663 ± 0.0006 GeV [2l|, while the prediction by 
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TABLE I: Total cross section for the hadronic production of -B^l 1 ^] an d B*[l Si] with varying 
m s , where nib = 4.9 GeV, ms s = m s + rrib, the gluon distribution function is taken from CTEQ6L, 
= Pt + m Bs an d a s is of leading order. 





TEVATRON( v / 5 = 1.96 TeV) 


LHC (VS = 14. TeV) 


m s (GeV) 


0.4 


0.5 


0.6 


0.7 


0.4 


0.5 


0.6 


0.7 


(T Bs {rib) 


48.25 


24.43 


14.33 


9.326 


512.6 


262.2 


155.1 


101.5 


a B * {nb) 


165.5 


82.18 


46.78 


29.30 


1739. 


871.8 


500.5 


316.2 




FIG. 2: B s differential distributions versus its Pt and rapidity y with varying m s . The gluon 
distribution function is taken from CTEQ6L, the characteristic energy scale is taken as Q 2 = 
p Tg + m Bs , and a s is of leading order. The solid line, dashed line, dotted line and dash-dot line 
stands for m s = 0.4, 0.5, 0.6 and 0.7 GeV respectively. The upper (lower) four lines corresponding 
to the distributions at LHC (TEVATRON) respectively. 



lattice QCD is about 5.37 GeV [22]. Thus with = 4.9 GeV and = + m s , the 
obtained is in the region of theoretical prediction as well as experimental measurement. 
In Table [H the total cross-section for the hadronic production of -Bsfl 1 ^] and -B*[l 3 5*i] at 
TEVATRON and LHC are computed. From TABHJ one may observe that the total cross 



2 This relation should be satisfied according to the gauge invariance of the hard subprocess. Further more, 
at the present, we treat B s as non-relativistic and weak binding state, then at LO the relative momentum 
between the constitute quarks can be ignored. 
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TABLE II: Total cross section for the hadronic production of -B^l 1 ^] and -B*[l 3 Si] with various 
mj,, where m s is fixed to be 0.5 GeV, the gluon distribution function is taken from CTEQ6L, the 
factorization energy scale is chosen = p\ + m% s and a s is of leading order. 





TEVATRON(^/5 = 1.96 TeV) 


LHC (VS = 14. TeV) 


m b (GeV) 


4.8 


4.9 


5.0 


4.8 


4.9 


5.0 


<TB s {nb) 


22.60 


24.43 


26.49 


280.9 


262.2 


245.7 


a B * (nb) 


88.86 


82.18 


76.09 


818.1 


871.8 


929.3 



section of ^[l 3 ^] is about 3 times bigger that of -B^l 1 ^], which roughly agree with the 
naive counting of spin states. m s affects the total cross section greatly, e.g. when m s increases 
by steps of 0.1 GeV, then the cross section of B s or B* decreases by about 80% — 100%. To 
show this point more clearly, we draw the B S -P F and rapidity Y distributions with m s = 0.4, 
0.5, 0.6 and 0.7 GeV respectively in FIGJ2J This implies that the present treatment of s- 
quark as heavy quark is reasonable but may be not too accurate, and a more accurate 
one, e.g. by including proper relativistic effects into the bound state, maybe improve the 
estimation, which is out of the range of the present letter. A similar calculation by varying 
mj, within its reasonable region but with fixed m s = 0.5 GeV, as shown by TABjTTl which 
shows that when mt, increases by steps of 0.1 GeV, the cross section of B s or B* changes 
slightly, which is around 10%. More precise values of and m s from potential model or 
lattice QCD can make the our estimations more reliable. In the following parts of the paper 
when examining the uncertainties from other factors, we shall always take the center values 
of m s = 0.5 GeV and = 4.9 GeV, for the quark masses. 

PDF is of non-perturbative nature, which can be obtained through global fitting of the 
experimental data. Here we take CTEQ6L and MRST2001lQ as typical examples to 
study the uncertainty caused by PDF. As shown in Eq.([T]) PDF can be factorized out at 
the energy scale [i 2 F with the help of pQCD factorization theorem. The factorization scale 
\j? F can be usually taken as the characteristic energy scale for the hard subprocess {Q 2 ), i.e. 
fip = Q 2 . For the present case, the gluon-gluon fusion subprocess is of three-body final state 
and contain heavy quarks, so there are ambiguities in choosing Q 2 and various choices of 
Q 2 would generate quite different results. Since such kind of ambiguity cannot be justified 
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TABLE III: Total cross-section for the hadronic production of i^fl 1 ^] and 2?*[l 3 5i] at TEVA- 
TRON and at LHC with LO running a s and the characteristic energy scale Type A: Q 2 = s/4; 
Type B: Q 2 = p| + m% and Type C: Q 2 = y 2 Th + mf . 





TEVATRON(\/S = 1.96 TeV) 


LRC(VS = 14. TeV) 




CTEQ6L 


CTEQ6L 


MRST2001 


CTEQ6L 


CTEQ6L 


CTEQ6L 


MRST2001L 


CTEQ6L 


Q 2 


Type A 


Type B 


TypeC 


Type A 


Type B 


Type C 


a Bs (nb) 


17.55 


24.43 


20.82 


25.01 


203.5 


262.3 


221.7 


264.1 


a B * {nb) 


58.37 


82.19 


69.83 


83.89 


675.8 


871.8 


734.2 


876.5 




FIG. 3: B s differential distributions versus its transverse momentum Pt and rapidity y for different 
LO PDFs, where the solid line and the dashed line are for CTEQ6L and MRST2001L respectively. 
The characteristic energy scale is taken as Q 2 = Pj> + m B s - The u PP er an d lower two lines 
corresponding to the distributions in LHC and TEVATRON accordingly. 

by the LO calculation itself, so we take it as the uncertainty of the LO calculation. In 
the following we choose three typical examples to study this kind of uncertainties: Type 
A: Q 2 = s/4, the CM. energy squared of the subprocess that is divided by 4; Type B: 
Q 2 — Pt + m B s > tli e transverse mass squared of the B s meson; and Type C: Q 2 = p^ b + ml, 
the transverse mass squared of the b quark. For comparison between TEVATRON and LHC 
and to pinpoint the uncertainties from PDFs, a s running and the choices of the characteristic 
energy scale Q 2 , we calculate the production cross sections according to two types of PDFs, 
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FIG. 4: B s differential distributions versus its transverse momentum px and rapidity y for typical 
choices of Q 2 , where the solid and the dashed lines are for Type A and Type B respectively. The 
gluon distribution is chosen as CTEQ6L and the running a s is in leading order. The upper and 
lower two lines corresponding to the distributions in LHC and TEVATRON accordingly. 

the strong coupling a s fixed by the corresponding PDFs and the characteristic Q 2 chosen as 
Type A, Type B and Type C. The obtained results are shown in TAB IIHI From TABLE |III| 
it is found that the difference caused by the two LO PDFs is small, which is ~ 15%. The 
choice of Q 2 as Type A and Type B cause changes is somewhat larger, i.e. 20% — 30%, 
while the choose of Type B and Type C leads to negligible changes to the cross section 
(less than 1%). The total cross section of the B(*> production at LHC are at least one 
order larger in magnitude than that at TEVATRON. This is mainly due to the fact that 
LHC (v 7 ^ = 14. TeV) has much higher collide energy than TEVATRON (v 7 ^ = 1.96 TeV), 
so the lowest boundary of the gluon momentum fractions X{ (i — 1,2) at LHC are much 
smaller than that at TEVATRON and then there are more interacting gluons that have a 
CM. energy above the threshold for the subprocess, in the collision hadrons at LHC than 
at TEVATRON. This can be shown more clearly by the Pt and y differential cross sections. 
More explicitly, we draw the curves for pseudo-scalar meson B s in FIGs.( )3|4l) . FIG J3] shows 
that the differential distributions for the two PDFs CTEQ6L and MRST2001L and FIGH 
shows that the differential distributions for the two Q 2 Type A and Type B. In regions of 
comparatively small pt and \y\, the distributions of MSRT2001L are smaller than that of 
CTEQ6L. From the figure, we also see that the pt distributions in TEVATRON are steeper 
than those in LHC. From TABLE IIHI we know that changes in the cross section caused by 
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TABLE IV: Dependence of R = ( sxbxxbqk ) on P Tcut for BJ^So] and B*[l 3 Si]. 

\ &LHC / 







5* 


Pr c «t(GeV) 


5 


20 


35 


50 





5 


20 


35 


50 


i?(xl0~ 2 ) 


9.32 7.71 


3.00 


1.67 


0.70 


9.43 


7.80 


3.12 


1.85 


0.87 




FIG. 5: B s differential distributions versus its y with various PTcut m TEVATRON (left diagram) 
and in LHC (right diagram). Solid line corresponds to the full production without PTcut', dashed 
line to PTcut = 5.0 GeV; dot line to prcut = 20.0 GeV; the dash-dot line to prcut = 35.0 GeV; the 
short dash line to PTcut = 50.0 GeV . 

Type B and Type C is quite small (less than 1%), and the curves of the production obtained 
by Type B and Type C are almost overlap, so we do not draw the curves for Type C. 

Experimentally, when the produced B s and B* mesons with a small Pt or a large rapidity 
y are too close to the collision beam, they cannot be measured, so only 'detectable' events 
should be taken into account, i.e. events with proper kinematic cuts on Pr and y should 
be properly set in the estimates. As a comparison, we define a ratio R = ( crTg J L ^ c f 0iV ) p 
to show how Pr C ut affects the integrated cross sections at TEVATRON and LHC, and the 
results is shown in TAB U VI It can be found that without Prcut, the integrated cross section 
at LHC is about one order higher than that at TEVATRON, and the value of R decreases 
greatly with the increment of Prcut, at about Prcut — 45 GeV, R ~ 1. 

Next, as an explicit example to show how the different cuts affect the production, we study 
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TABLE V: Values of Rp Tcut for the hadronic production of pseudo-scalar B s meson in TEVATRON 
and LHC. 



PTcut 


0.0 GeV 


5 GeV 


20 GeV 


35 GeV 


50 GeV 


Vcut 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


R PTcut (TEVATRON) 


0.47 


0.66 


0.80 


0.48 


0.67 


0.81 


0.59 


0.79 


0.92 


0.68 


0.86 


0.95 


0.69 


0.92 


0.99 


R PTcut (LHC) 


0.32 


0.47 


0.60 


0.33 


0.48 


0.61 


0.40 


0.57 


0.71 


0.45 


0.62 


0.77 


0.56 


0.73 


0.83 



the distributions of Pt and y for B s . For the present purpose, we take CTEQ6L for PDF, 
LO running a s and Type B energy scale to carry out the study. The correlations between 
Pt and y are interesting, so we plot the ^-distributions with various Pt-cvXs over a wide 
range Pxcut '■ 5.0 ~ 35 GeV in FIGJ5J FIGJ5] shows that the dependence of the differential 
distributions on rapidity y with different PTcut at LHC exhibits a broader profile than that 
at TEVATRON. The p^-distributions of the production vary with y cut mainly due to the fact 
that as pt increases, the dependence of the distribution on y becomes smaller as the value 
of ycut becomes less important. To analyze the quantitative difference of the differential 
distributions with regard to PTcut and y cut , we introduce a ratio for the integrated hadronic 
cross sections, Rp Tcut = (^^J , where a ycut and o"o are the hadronic cross section with 
and without y cut respectively. The ratio R PTcut varies with PTcut and y cut , and its values are 
given in TAB|V] TAB JVl shows that for a fixed y cut , the value of Rp Tcut becomes larger with 
increasing PTcut- It is understandable that the differential distributions versus the rapidity 
y decrease with the increment of Pt, so the contributions to the hadronic cross section 
surviving after the cut, i.e. (\y\ < y cu t), increase with the increment of PTcut- 

To summarize: We have presented quantitative studies on the uncertainties in estimates of 
the B^ meson hadronic production within the FFN scheme. The investigated quantitatively 
uncertainties involve the PDF, the values of rrib and m s , and the characteristic energy scale 
Q 2 of the process and etc.. It is found that when m s increases by steps of 0.1 GeV, the 
integrated cross section of B^ decreases by about 80% — 100%, when increases by steps of 
0.1 GeV, it is about 10%. While the uncertainties caused by the parton distribution function 
and the factorization scale varies within the region of | to |. We have also shown the 
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differences between LHC and TEVATRON for various observable with reasonable kinematic 
cuts, such as the cuts on the B^ meson transverse momentum Pxcut and rapidity y cut . 
Our results show that the experimental studies of the Bj> meson at the two colliders are 
complimentary and stimulative. Concerning the prospects of B s production at Fermilab 
TEVATRON and at CERN LHC, the obtained results may be as useful references for these 
experiments. Since LHC has much higher luminosity and higher collision energy than that of 
TEVATRON, it seems that the particularly interesting topics on B s may be more accessible 
and fruitful at LHC than that at TEVATRON. Further more, it is reasonable to assume 



that, similar to the hadronic production of J/^, B c and H cc |25j, [26|, l27(], the 'heavy quark 
mechanisms' via the sub-processes g + s — > 

SW + ... and g + b -> flW + ... may be as 
important as the gluon-gluon fusion mechanism, which should be treated on the equal footing 
in comparison to that of the gluon-gluon fusion mechanism. However to be consistent 
theoretically and to deal with the possible double counting from all these mechanisms, 



one should work in the general-mass variable-flavor-number (GM-VFN) scheme 28 



23, 



301 ] in stead of the FFN scheme. A detailed discussion on the GM-VFN scheme, and a 



comparison of B^ production within these two schemes is in preparation and shall be 
presented elsewhere. 
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